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Abstract 

The lac operon of Escherichia coli can exhibit bistability. Early studies showed that bistability occurs dur- 
ing growth on TMG/succinate and lactose + glucose, but not during growth on lactose. More recently, 
studies with lacGFP-transfected cells show bistability during growth on TMG/succinate, but not during 
growth on lactose and lactose + glucose. In the literature, these results are invariably attributed to varia- 
tions in the destabilizing effect of the positive feedback generated by induction. Specifically, during growth 
on TMG/succinate, lac induction generates strong positive feedback because the permease stimulates the 
accumulation of intracellular TMG, which, in turn, promotes the synthesis of even more permease. This 
positive feedback is attenuated during growth on lactose because hydrolysis of intracellular lactose by (3- 
galactosidase suppresses the stimulatory effect of the permease. It is attenuated even more during growth 
on lactose + glucose because glucose inhibits the uptake of lactose. But it is clear that the stabilizing effect 
of dilution also changes dramatically as a function of the medium composition. For instance, during growth 
on TMG/succinate, the dilution rate of lac permease is proportional to its activity, e, because the specific 
growth rate is independent of e (it is completely determined by the concentration of succinate). However, 
during growth on lactose, the dilution rate of the permease is proportional to because the specific growth 
rate is proportional to the specific lactose uptake rate, which in turn, proportional to e. We show that: 
(a) This dependence on creates such a strong stabilizing effect that bistability is virtually impossible 
during growth on lactose, even in the face of the intense positive feedback generated by induction, (b) This 
stabilizing effect is weakened during growth on lactose + glucose because the specific growth rate on glucose 
is independent of e, so that the dilution rate once again contains a term that is proportional to e. These 
results imply that the lac operon is much more prone to bistability if the medium contain carbon sources 
that cannot be metabolized by the lac enzymes, e.g., succinate during growth on TMG/succinate and glucose 
during growth on lactose + glucose. We discuss the experimental data in the light of these results. 
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1 Introduction 

The lac operon has been a topic of considerable interest since the late 1940's. This interest was 
stimulated by the hope that insights into the mechanism of lac induction would shed light on the 
central problem of d evelopm e nt, namely, the mechani sm by which genetically identical cells acquire 
distinct phenotypes ( Monodl . 1947 : Spiegelman . 19481 ). 



Many of the early studies were concerned with the kinetics of enzyme induction. Initial attempts 
to measure the kinetics were hindered by the fact that lactose, the substrate that stimulates the 
induction of the lac operon, promotes not only the synthesis of the lac enzymes, but also their dilu- 
tion by growth. Under these conditions, it is impossible to separate the kinetics of enzyme synthesis 
from the masking effects of dilution. This obstacle was overcome by the discovery of gratuitous in- 
ducers, such as methyl galactoside (MG) and thiomethyl galactoside (TMG). Enzyme synthesis and 
dilution could be uncoupled by exposing the cells to a medium containing a gratuituous inducer and 
non-galactosidic carbon sources, such as glucose or/and succinate. The gratuitous inducer promoted 
enzyme synthesis, but not growth, and the non-galactosidic carbon sources supported growth, but 
not enzyme synthesis. 

Although gratuitous inducers enabled enzyme synthesis and dilution to be uncoupled, an impor- 
tant question remained. Specifically, it was not known whether gratuitous inducers provoked the 
same enzyme synthesis rate in every cell of a cultu re. Initial experiments suggested that this was 
indeed the case. Benzer showed that ( Benzer . 19531 . Figs. 6 and 7): 



1. If non-induced cells of E. coli B (pregrown on lactate) were exposed to 1 g/L of lactose, 
only a small fraction of the cells synthesized /3-galactosidase initially. This fraction increased 
progressively until the culture became homogeneous eventually. 

2. In sharp contrast, if non-induced cells were exposed to 2 g/L of MG, all the cells started syn- 
thesizing /?-galactosidase immediately and at near-maximal rates, i.e., the population became 
homogeneous almost instantly. 

However, it was shown later that the population became homogeneous instantly only because the 
concentration of the gratuitous inducer was high. At low concentrations of the gratuitious inducer: 

1. The population remained heterogeneous for a significant period of time. Furthermore, the 
smaller the concentration of the gr atuitous inducer, the longe r the time required for the 



population to become homogeneous (ICohn and Horibatal . Il959al . Table 1 



2. The enzyme synthesis rate was not uniquely determined by the composition of the medium 
(Fig. dK). If TMG and glucose were added simultaneously to a culture of E. coli ML30 
growing on succinate, there was almost no synthesis of /?-galactosidase. However, if TMG 
was added 15 mins before the addition of glucose, /3-galactosidase was synthesized for up to 
130 generations. Thus, enzyme synthesis is bistable: Pre-induced cells remain induced, and 
non-induced cells remain non-induced. 

The existence of bistability and heterogeneity depended crucially upon the existen ce of lac permease 



(LacY ). Both phenomena disappeared in lacY (cryptic), but not lacZ , mutants (|Cohn and Horibata 
19593). 



Subsequent experiments showed that bistability also occurred during growth of E. coli K12 3.000 
on lactose + glucose (Fig. [T]b). If lactose was added to a culture growing on glucose, there was no 
/3-galactosidase synthesis. However, if the culture was exposed to 10~^ M IPTG before the addition 
of lactose, /3-galactosidase synthesis persisted for several generations. 
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Fig. 1: 



Bistability duri ng growth of E. coli on TMG (jCohn and Horibatal . Il959bl . Fig. 4) and lac- 
tose + glucose (jLoomis and Magasanikl . 119671 . Fig. 3). (a) If glucose and TMG are added 
simultaneously to a culture growing on succinate, there is no /3-galactosidase synthesis (♦). If 
glucose is added to the culture 20 mins after the addition of TMG, the enzyme is synthesized 
(■) at a rate that is 50% of the rate observed when only TMG is added to the culture (A), 
(b) If lactose is added to a culture growing on glucose, there is no /3-galactosidase synthesis 
(♦). If lactose and glucose are added to the culture after it has been exposed to IPTG for 
15 mins, the enzyme synthesis rate (■) increases within a few generation to ~40% of the 
rate in a culture exposed to IPTG only (A). 



m 



These intrica t e dyn amics attracted significant attention among some theoreticians (reviewed 



Laurent et all . l2005l l. In particular, Babloyantz & Sanglier formulated a model of growth on 



TMG/succinate which took due ac count of enzyme synthesis by th e Jacob-Monod mechanism, and 
enzyme depletion by degradation (iBabloyantz and Sang ii^ . I1972I ). They showed that the model 
yielded the bistability observed in experiments. Chung and Stephanopoulos formulated a simi- 
lar model, the main differences being that repressor-operator and repressor-inducer binding were 
assu nied to be in quasi-equilibrium, an d the enzyme was depleted by both degradation and dilu- 



tion ( Chung and Stephanopoulos . 19961 ). This model is given by the equations 



dx 
'dt 
de 
di 



Ts-r^ - rgX, Ts = VgC 



l + Kix 



2J2 



1 + Q + &2' ^ 



T ^ - — - /c^ 6 



(1) 

(2) 



where x and s denote the intracellular and extracellular TMG concentrations, respectively; e denotes 
the lac permease activity; Vg is the specific growth rate on the non-galactosidic carbon source; rs,r~ 
denote the specific rates of TMG uptake and expulsion, respectively; and rf,r~ denote the specific 
rates of permease synthesis and degradation, respectively. The expression for r+ is based on the 
molecular model formulated by Yagil & Yagil, which assumes that the l ac operon conta i ns on e 
operator, and the lac repressor contains identical two inducer-binding sites ( Yagil and Yagij 197ll ). 
The parameter, Kx, is the association constant for the repressor-inducer binding; and a is jointly 
proportional to the intracellular repressor level and the association constant for repressor-operator 
binding. Evidently, a is a measure of the repression, defined as the ratio, r+l^^^/ r+|^^Q. 

Although the experiments done by Cohn and coworkers provided clear evidence of bistability 
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Fig. 2: Dynamics of the lac operon (from lOzbudak et"al1 . 120041 1 . (a) Bistability in wild-type cells. 



The (normalized) green fluorescence provides a measure of the steady state activity of the 
lac operon. The upper (resp., lower) panel shows the green fluorescence observed when 
an induced (resp., non-induced) inoculum of Escherichia coli is grown exponentially on a 
mixture of succinate and various concentrations of extracellular TMG. Bistability occurs at 
TMG concentrations between 3 and 30 ^M: The green fluorescence is ~100 if the inoculum is 
fully induced, and ~0.5 if the inoculum is non-induced, (b) Monostability in low-repression 
cells exposed to 1 mM glucose and various concentrations of extracellular TMG. 



during growth on TMG/succinate, they did not investigate the enzyme levels at a wide variety of 
conditions. Re cently, Qzbudak et al measured the steady state enzyme levels at various concentra- 



tions of TMG JOzbudak et a,l.l . l2004l l. To this end, they inserted into the chromosome of Escherichia 



coli MG 1655 a reporter lac operon, i.e., an operon under the control of the lac promoter, which 
codes for the green fluorescent protein (GFP) instead of the lac enzymes. They then exposed non- 
induced and induced cells to a fixed concentration of succinate, and various concentrations of TMG. 
It was observed that: 



1. When the cells are grown in the presence of succinate and various concentrations of TMG, 
they exhibit bistability (Fig. [2^). This bistability persists even if glucose is added to the 
mixture of succinate and TMG. 



2. The bistability disappears if the concentration of the lac repressor is reduced 
transfecting the cells with the lac operator (Fig. [2b). 



-40-fold by 



They also showed that these two observations were mirrored by the bifurcation diagram for the 
Chung-Stephanopoulos model. 

The Yagil & Yagil model of lac induction is not consistent with the structure of the lac operon 
and repressor. The lac operon contains two auxiliary operators, O2 and O.s, in addit ion to the main 
operator, Oi, and the lac repressor contains four inducer-binding sites ( Lewis . 20051 ). Furthermore, 
these structural f eatures play a cruc i al rol e in the formati on of DNA loops, th e key determinants 
of lac repression JOehler et al.l . Il990l . flOoi ) and induction (jOehler et al.l . l200d l. Molecular models 
taking due account of the 3 operators and 4 inducer-binding sites yield the lac induction rate 



'l + a/{l + K^xy + a/{l + K^x) 



4 ' 



(3) 



where Kx is the association constant for repressor-inducer binding, and a, a are rela ted to the lac re- 
pression stemming from repressor-operator binding and DNA looping, respectively (jKuhlman et al.l . 
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Fig. 3: The dynamics of the inducer a r e diffe rent during growth on (a) TMG/succinate and (b) lac- 
tose (adapted from ISavagead . l2nnih . Here, 5 denotes extracellular TMG or lactose; X 
denotes intracellular TMG or lactose; M denotes lac mRNA; and Ep,Eb denote lac perme- 
ase and /9-galactosidase, respectively, (a) Intracellular TMG is discharged from the cell by 
inducer expulsion, a process that is independent of lac enzymes, (b) Intracellular lactose is 
metabolised by /3-galactosidase {Eb). 



20071 : iNarang) . 120071 : ISantillan et al.L 120071 ). In wild-type lac, the repression, 1 -|-a-|-a, is 1300, and 
the bulk of this repression is due to DN A looping (a !ii 20, a 12 50) mediated by the interaction 



of repressor- bound Oi with O2 and O3 (lOehler et al.l . Il990l . 1 19941 ). The first goal of this work is 
to determine if the dynamics of the Chung-Stephanopoulos model are significantly altered by these 
more realistic kinetics. To this end, we consider the modified Chung-Stephanopoulos model in which 
the induction rate is replaced by eq. JS]). We show that the dynamics of this modified model are in 
quantitative agreement with the data. 



Ozbudak et al also studied the growth on lactose and lactose + glucose (iQzbudak et al.l . 12004 



p. 2 of Supplement). They found that when non-induced cells (pregrown on succinate) are exposed 
to various concentrations of lactose and lactose + glucose, the green fluorescence of the cells has a 
unimodal distribution after 4 hours of growth. They did not report any experiments with induced 
cells. However, the data for TMG/succinate shows that the green fluorescence of non-induced cells 
has a bimodal distribution near the upper limit of the bistable region (corresponding to extracellular 
TMG levels of 15-30 ^M in Fig. [2^). The absence of such a bimodal distribution led them to conclude 
that bistability does not occur during growth on lactose and lactose -|- glucose. This is consistent 
with the data obtained by Benzer, but contradicts the data shown in Fig. [lb. The second goal of 
this work is to seek an explanation for the absence of bistability during growth on lactose, and the 
conflicting results for lactose + glucose. 

These experimental results have spurred the development of several mathematical models, most 
of which are concerned with the disappearance of bistability during growth on lactose. Thus far, 
two mechanisms have been proposed. 

The first mechanism proposes that during growth on lactose, the induced cells outgrow the non- 
induced cells. It seems unlikely that this mechanism, by itself, can explain the data. To see this, 
suppose that the non-induced inoculum used in the experiments contains 10% induced cells, which 
double every hour. As a worst-case scenario, assume that the non-induced cells do not grow at all. 
Then, after 4 h, almost 40% of the population is still non-induced, which is far from the unimodal 
distribution observed in the experiments. 

The second mechanism proposes that bistability does not eve n exist because positive feed- 
back is suppressed during growth on lac t ose or lactose + gluc ose (jvan Hoek and Hogeweg) . l2006l : 
Mahaffy and Savevl . I1999I : ISantillan et al.l . l2007l : IbavageatJ . l200lh . More precisely, during growth on 
TMG/succinate, bistability is feasible because of strong positive feedback: The permease stimulates 
the accumulation of intracellular TMG, which in turn promotes the synthesis of even more permease 
(Fig. [3tt). The destabilizing effect of this positive feedback produces bistability by overcoming the 
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stabilizing effect of dilution. During growth on lactose, the positive feedback is suppressed because 
hydrolysis of lactose by /3-galactosidase attenuates the stimulatory effect of the permease (Fig. [3b)- 
It is attenuated even more during growth on lactose + glucose because in the pres ence of glucose , 



enzym e IIA^'*^ is dephosphorylated, and inhibits the permease by binding to it (jSantillan et al 



20071). 

In all these models, the differences in the dynamics of growth on TMG/succinate, lactose, and 
lactose + glucose are attributed entirely to changes in the destabilizing effect of positive feedback 
generated by induction. Here, we show that the stabilizing effect of dilution also changes dramat- 
ically with the medium composition, and this has equally profound effects on the dynamics of the 
lac operon. Specifically: 

1. The stabilizing effect of dilution is much stronger during growth on lactose (as opposed to 
growth on TMG/succinate). Indeed, during growth on TMG/succinate, the dilution rate of 
the lac enzymes is proportional to their level, e, because the specific growth rate does not 
depend on the activity of these enzymes — it is completely determined by the concentration 
of succinate. However, during growth on lactose, the dilution rate is proportional to be- 
cause the specific growth rate is proportional to the lactose uptake rate, which, in turn, is 
proportional to the activity of lac permease. We show that this stronger stabilizing effect of 
dilution suppresses bistability on lactose even in the presence of the intense positive feedback. 

2. The enhanced stabilizing effect of dilution is attenuated once again during growth on lactose + 
glucose. This is because the specific growth rate on glucose, a non-galactosidic carbon source 
like succinate, is independent of e. Thus, in the presence of glucose, the dilution rate once 
again contains a term that is proportional to e, and the dynamics become similar to those on 
TMG/succinate, i.e., bistability is feasible, provided the positive feedback is sufficiently large. 

These results imply that the lac operon is much more susceptible to bistability in the presence of 
non-galactosidic carbon sources, since they serve to suppress the stabilizing effect of dilution. 

2 The model 

Fig. m shows the kinetic scheme of the model. Here, Si, S2 denote the exogenous lactose (or TMG) 
and glucose, respectively; Ei,E2 denote the transport enzymes for lactose (or TMG) and glucose, 
respectively; Xi,X2 denote internalized lactose (or TMG) and glucose, respectively; G denotes the 
GFP synthesized by the reporter lac operon; and C~ denotes all intracellular components except 
Ei, Xi, and G (thus, it includes precursors, free amino acids, and macromolecules) . We assume 
that: 

1. The concentrations of the intracellular components, denoted Cj, Xj, g, and c~, are based on 
the dry weight of the cells (g per g dry weight of cells, i.e., g gdw~^). The concentrations of 
the exogenous substrate and cells, denoted and c, are based on the volume of the reactor 
(g/L and gdw/L, respectively). The rates of all the processes are based on the dry weight of 
the cells (g gdw~^ h~^). We shall use the term specific rate to emphasize this point. 
The choice of these units implies that if the concentration of any intracellular component, Z , 
IS z g gdw~^, then the evolution of z in batch cultures is given by 

dz _|_ _ f 1 dc 
It - '^z - [^-^ 

where rf and r~ denote the specific rates of synthesis and degradation of Z in g gdw"^ h^^. 
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Fig. 4: Kinetic scheme of the model. 



2. The specific uptake rate of Si, denoted r<j^j, follows the modified Michaelis-Menten kinetics, 

fs,! = ^s,i(^iSi / {Kg^i + Sj). 

3. In the case of glucose or lactose, part of the internalized substrate, denoted Xj, is expelled 
into the environment. The remainder is converted to C~, and oxidized to CO2 for generating 
energy. 

(a) The specific rate of expulsion of Xi follows first-order kinetics, i.e., r~- = k~-Xi. 

(b) The conversion of Xi to C~ and CO2 also follows first-order kinetics, i.e., rj"^ = k^iXi. 

(c) The specific rate of synthesis of C~ from Xi is Yir^^, where Yi is a constant (which will 
turn out later to be effectively equal to the yield of biomass on Si). 

4. Internalized TMG is completely expelled — it does not support biosynthesis or respiration. 

5. The internalized substrates induce the synthesis of the enzymes and GFP. 

(a) The specific synthesis rate of the lactose enzymes, Ei, follows the kinetics 



reA = Ve 



1 



e,l 



1 + ai/ (1 + K^^iXiY + di/ (1 + 



(4) 



where K^^i is the association constant for repressor-inducer binding, and ai,ai charac- 
terize the transcriptional repression in the absence of inducer due to repressor-operator 
binding and DNA loop ing, respective lv- Both ai and di are proportional to the intra- 
cellular repressor level (jNarang) . 120071 ) . 



In wild-type cells, ai = 20, di = 1250 (lOehler et al.l . ll990, Fig. 2). In ce l ls tra nsfected 



with the lac operator, the repressor levels decrease 43-fold 
ai 0, di « 30. 



Ozbudak et"al] . l2004l ): hence 
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(b) The specific synthesis rate of GFP follows the same kinetics as the lac operon, i.e., 

1 



rGFP = Vgfp 



1 + ac/ (1 + Kc^^ixiY + qg/ (1 + Kccixi] 



We shall assume that the promoters of the reporter and native lac operons are identical, 
so that Vgfp = Ve,! and = ai = 20. However, ac < Oii because the reporter lac 
operon lacks the auxiliary operator, O2, which precludes the formation of DNA loops 
due to interaction between Oi and O9,. Ozbudak et al found the repression of the lac 



reporter, l + ac + ac to be 170 (lOzbudak et al.l . 120041 ) . which implies that ac ~ 150. If 



the cells are transfected with the lac operator, ac ~ 0, ac = 4. 

(c) The induction of the glucose enzvmes , E2, is thought to occur by a mechanism similar to 
the one that induces the lac operon (|Plumbridgel . l2003l l. Specifically, in the absence of 



glucose, transcription of the ptsG operon is blocked because the repressor (Mlc) is bound 
to the operator. In the presence of glucose, the enzyme IF^*^ sequesters Mlc from the 
operator by an unknown mechanism, thus liberating the operon for transcription. We 
assume that the specific synthesis rate of E2 has the form 



^e,2 = VeX 



where Kx^2 and 02 are phenomenological parameters, i.e., they cannot be expressed in 
terms of parameters characterizing the molecular interactions. 

The data shows that the PTS level in cells growing exponentially on glucose is roughly 5 
times the level observed in cells growing exponentially on glycerol (Fig. [5]). This implies 
that a2 ~ 4. 

(d) The synthesis of the enzymes and GFP occurs at the expense of the biosynthetic con- 
stituents, C~ . 

(e) Enzyme and GFP degradation are negligibly small. 

(f) Non-specific diffusion of the substrates into the cell is negligibly small. This is valid 
for lactose and glucose at the concentrations typically used in the experiments. It is 
valid for gratuitous induce rs, such as TMG, only if the extracellular concentration is 
<50 //M |Herzenbergl . ll959l . Fig. 4). 



In what follows, we begin by deriving the equations for growth on TMG plus non-galactosidic 
substrates. We then derive the equations for growth on lactose + glucose, from which the equations 
for growth on lactose are obtained by letting the concentration of glucose be zero. 
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Fig. 5: The phosphotransferase system of enzymes is inducible ( Bettenbrock et al. . 20061 ). During 
batch growth of E. coli on glycerol, the PTS activity (x) is ~ 0.01. The concentrations of 
cells (♦) and glucose (a) correspond to the growth on glucose + glycerol. During the first 
phase of diauxic growth on glucose + glycerol, the PTS activity (■) increases to ~0.05. 



2.1 Growth in the presence of TMG and non-galactosidic carbon sources 

During growth in the presence of TMG and non-galactosidic carbon sources, the mass balances yield 

si 



dsi 
If 
dxi 
IT 
dei 
If 

dg_ 
dt 

dc 
di 



- Vs,iei 

Ve,l 
Ve,l 



Ks,i + si 

Sl 



KsA + Sl 



1 



1 + Qi/ (1 + K^,ixif + ail (1 + K^,ixif 



1 + ai/ (1 + K^,ixiy + acl (1 + K^,ixi] 



Idc 

-cdft^''' 
1 dc 

-cdft^'' 



r„c 



where r^, the exponential growth rate on the non-galactosidic carbon sources, is a fixed parameter 
(independent of the model variables, si,xi,ei, g,c). It is completely determined by the concentration 
of the non-galactosidic carbon sources(s). 

It turns out that the dynamics of the experiments can be described by only 2 differential equa- 
tions. Indeed, since the cell density remains vanishingly small throughout the experiment, there is 
almost no depletion of extracellul ar TMG, and si remains essentially equal to its initial value, si^o- 



Moreover, since A;^ ~ 1 min ^ (IKeped . Il960l . Table 1) and ~ 0.5 hr ^, xi attains quasisteady 



state on a time scale of minutes. It follows that the dynamics of Ei after the first few minutes are 
well- approximated by the equations 



dei 
If 



1 



e,l 



" 1 + ai/ (1 + K^^ixif + ail (1 + K^^ixif 



da _ 

— — Ve.l 



dt 

Xl 



' 1 + ai/ {1 + K^^iXiY + ac/ (l + K^^ixi) 



sA 



-eicJi, cJi 



^1,0 



(5) 
(6) 

(7) 
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where (l7| shows the quasisteady state concentration of intracellular TMG. 

Eqs. ([5|)-(l6]) imply that the maximum steady state level of Ei and G is Ve^i/rg. It is therefore 
natural to define the dimensionless variables 



_ei ^ 5^ 



ei = TP — J—' 7 = — — Xi = K^^ixi, T = Tgt. 



This yields the dimensionless equations 

dei 1 



dT i + ail {l + xif + ai/ {l + xif 

^ = I ^ (9) 

dr i + a,/{l + xif + aG/{l + Xif 

XI = (10) 

where 

C>1 = C>l,mO"l, = — 1 ■ (11) 

It follows from (fTO]) - (fTT]l that Si^m is a measure of the quasisteady state intracellular TMG con- 
centration at saturating levels of extracellular TMG (measured in units of K^\, the dissociation 
constant for repressor-inducer binding). We can also view 5i^ra as a measure of the strength of the 
positive feedback: Given any di, the strength of the positive feedback, dve^i/dei, is an increasing 
function of 6i^m- The parameter, is a measure of the intracellular TMG concentration at any 
given extracellular TMG level. For a given cell type and non-galactosidic carbon source, 5i is pro- 
portional to (Ti, and hence, can be treated as a surrogate for the extracellular TMG concentration. 

2.2 Growth on lactose or lactose + glucose 

During growth in the presence of lactose and glucose, the mass balances yield 

dxi Sj + fldc\ 

dt Ks^i + Si ""'^ \c dt J 

— V '^'^^ e (li) 

1 + ai/ (1 + ir.,ixi)' + ai/ (1 + K^^^xif ~ \c1t) 

dt '•'l + ai/(l + i^,,iXi)VaG/(l + i^.MXi)' \cdtj^' 

^^1/ I ( l^A Hfi^ 

dc~ v^/, ^,-1- \ fldc\ 



dt E - - rCFP - re,2 " J) (17) 

i=l 
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It is shown in Appendix [A] that under the experimental conditions, the dynamics of the enzymes 
and GFP are well- approximated by the equations 

-IT = ^e,i- --; — -; -4 - {(piYiVs^iaiei + </'2>2^s,2'72e2) ei, (18) 

-j: = --; ^ — -; ;4 - (^lyil^.i^iei + '/'2>2 ^5,20-262) g, (19) 

QjG'? 1 

-JT = Ve,2— TTm? T - + <p2Y2Vs,2(T2e2) 62, (20) 

dt 1 + 02/ (1 + Kx,2X2) 

•^i — TTjl _ SiCTi, (21) 



x,t ' x,t 



where 



Ks,i + Si^o ' * k^.-\-k^^ 



The parameter, (pi, is the fraction of substrate intake that is channeled into growth and respiration 
(the remainder is expelled into the medium). The second term in eqs. (fT8]) - ([20]) represents the 
dilution rate of Ei, G, and E2, respectively!! 

It follows from (fT8]) - ([20]) that during single-substrate growth on Si, the steady state activity of 
Ei is at most 



y (piYiVs^iai' 
and the maximum specific growth rate is at most 



^i^s,i\l ^ ^I'l 0-j = ^J 4>iYiVs,iVe,i(Ti. 



Thus, we are led to define the dimensionless variables 

_ e.i _ 9 



which yield the dimensionless equations 

dei 1 



,7 = /ta //i t^t^ ^« = K^^iXi, T = t^J (t)iYiVs,lVe,l(Tl, 



dT i + ai/ {l + xif + ai/ (l + xi)^ 
drf_ 1 

dT~ 1 + axl (1 + xi)^ + olgI (1 + Xi)^ 



(ei+ae2)ei, (22) 



(ei+ae2)7, (23) 



dt2 1 / \ 

dT l + a2/(l + X2) 

Xi = kei,i = \,2 (25) 
with dimensionless parameters 



^ ^ a/0212 V's,2Ve,2tT2 ^^g) 
y^0lYlFs,lVe,iCri ' 

<Ji = Oi,my Gi, Oira = — 7^^ \\ — TTJ~ ■ l^'J 



^ These equations are formally similar to the model considered in lNarang and Pilvuginl . l2007l . the main difference 
being that the induction kinetics were assumed to follow Yagil & Yagil kinetics. 
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Here, a, is a measure of the specific growth rate on S2 relative to that on Si, and Si^m is a measure 
of the quasisteady state concentration of at saturating concentrations of Si (or equivalently, the 
strength of the positive feedback generated by induction of Ei). 

Unhke TMG, lactose is rapidly metabolized to support growth and respiration. It follows that 
the ability of the cells to accumulate the intracellular substrate, (and hence, the strength of the 
positive feedback) is smaller during growth on lactose, i.e., 5i^m < ^i,m- Indeed, (fTT]l and ([27l) imply 
that 



Since r^, the specific growth rate in the experiments with TMG, is comparable to y^^i^iT^^i^^^, a 
measure of the maximum specific growth rate on lactose, we have 

which is less than 1. 

In the particular case of growth on lactose, (T2 = a = 0, and the above equations become 



dei 1 2 



dr i + a^/(i + ^^)2 + ai/(l + xi) 
dj 1 



ei, (28) 



ei7, (29) 



dT i + ai/il + xif + ac/il + Xif 
XI = hei. (30) 

Note that ((28]) is formally similar to jS]), the only difference being that the dilution rate is propor- 
tional to e\ rather than ei. This reflects the fact that during growth on lactose, the specific growth 
rate is proportional to the activity of Ei. 

3 Results and Discussion 

We note at the outset that the steady state GFP level is completely determined by the steady state 
activity of the lac enzymes. More precisely, eqs. ©-([Q]) and ([22]) - ([23]) imply that 

2 ^ l + ai/(l + Xi)' + «i/(l + Xi)' .3^. 

ei l + ai/(l + xi)' + aG/(l + Xi)"' 
where xi is given by (fTOl) or ([25]) . In what follows, we shall focus on the variation of the steady 
state enzyme activity, ei, with the extracellular TMG (5i) or lactose (5i) concentration. Given this 
relation, the steady state GFP level completely determined by ([3T]l . 

If the repression characteristics of the lac reporter were identical to those of native lac (ai = ac, 
ai = ac), 7 would be identical to ei. However, since ac = 150 is significantly smaller than 
ai = 1250, the ratio, 7/ei, is a decreasing function of xi- Now the inducer levels, xi, are vanishingly 
small in non-induced cells, and very large in induced cells. Hence, ([3TI1 implies that 

Tnon-induced ^ Tinduced ^, ^ ^ (Tinduced/Tnon-induced) 1 



,non-induced ^1, induced V^l, induced /ei,non-induced) 8' 

i.e., ratio of the GFP levels in induced and non-induced cells is significantly smaller than the 
corresponding ratio of the enzyme activities. We shall appeal to this fact later. 

In what follows, we consider the growth on TMG, lactose, and lactose -|- glucose. In the first 
two cases, we also study the dynamics in the absence of DNA looping {ai = 0). We consider this 
biologically unrealistic scenario because it yields useful intuitive insights. 
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Fig. 6: Variation of the steady state enzyme activity (ei) with extracellular TMG level (^i) at high 
and low repression levels, (a) At high repression levels (ai = 1200), there is a range of 
extracellular TMG concentrations at which the enzyme activity is bistable, (b) At low 
repression levels (ai = 25), there is a unique enzyme activity at every extracellular TMG 
level. 



3.1 Growth on TMG and non-galactosidic carbon sources 
3.1.1 No DNA looping 

In this case, the enzyme dynamics are given by ([8ll- (fTQ]l with qi = 0, and the steady states satisfy 
the equation 

/(^i) ^ , ,2 - = 0- (32) 

1 + ai/ {l + 5iei) 

Since the induction rate lies between 1/(1 + ai) and 1, so does the steady state enzyme activity. 
Eq. ([32]l captures the steady state data shown in Fig. [2l Indeed, ([32]l implies that 



<3i(ei 



ei (1 - ei) ei 

Given any ai, the parametric curve, (5i(ei),ei), 1/(1 + ai) < ei < 1, yields the variation of the 
steady state enzyme activity with ^i, a surrogate for the extracellular TMG concentration. If the 
repression is large, the locus of steady states is hysteretic (Fig. [6|i). If the repression is small, there 
is a unique enzyme activity at every extracellular TMG level (Fig. [6b). 

The loss of bistability at low repression levels becomes more transparent if we plot the surface 
of steady states as a function of the parameters, ai,6i (Fig. [7h). The steady states in Fig. [6] were 
obtained by varying at fixed ai. These steady states are represented in Fig.[7h by the intersection 
of the plane, ai = constant, with the surface of steady states. Now, at high repression levels, the 
surface contains two folds (red and blue curves). Thus, the steady states of high-repression cells 
correspond to a curve with two folds. As the repression level decreases, the two folds approach each 
other until they merge and disappear. The steady states of low-repression cells therefore increase 
monotonically. 

Ozbudak et al showed that the transition from bistability to monostability occurs at a repression 
level lying between 5 and 50. We show below that the model captures this result. To this end, we 
begin by deriving the equations describing the folds of the steady state surface. Evidently, the folds 
satisfy ((321) because they lie on the steady state surface. Since the folds are singular points of the 
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Fig. 7: (a) The surface of steady state enzyme levels for growth on TMG. The surface contains two 
folds, which are represented by the red and blue curves, (b) The bifurcation diagram for 
growth on TMG, obtained by projecting the folds in (a) onto the ai^i-plane. Bistability 
occurs precisely if ai and lie in the region between the red and blue curves. The dashed 
line show the path along which 5i changes when a given cell type (fixed ai) is exposed to 
various extracellular TMG concentrations. 



surface, they also satisfy the equation 



2ai5i/(l + 5iei)' 



l + ai/(l + 5iei)' 



1 = 0. 



(33) 



It is shown in Appendix[B]that eqs. ([32]) - ([33]l define a curve in the ai^iei-space with the parametric 
representation 



ei(xi 



XI -1 



ai(xi 



(1+Xi) 
1 



2xi 



(34) 



2X1 " ' Xi - -L ' " ■ Xi 

As xi increases from 1+ to oo, this parametric representation traces the folds of the steady state 
surface. 

The bifurcation diagram for a system refers to the classification of its dynamics in parameter 
space. For the system at hand, it is obtained by projecting the folds onto the ai^i-plane (Fig. [Tb). 
Evidently, multiple steady states occur precisely when the parameters, ai and ^i, lie between the 
blue and red curves in Fig. [7b. These two curves meet at a cusp characterized by the conditions 



dai _^ (l + xi)'(xi-2) 
dxi ixi - if 

fi = 2^^(^^-f =0, 
^Xi (xi - 1)' 



0, 
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which imply that the cusp occurs at xi = 2, and its coordinates are ai = 27, 6i = 8. It follows 
that bistability is feasible only if ai > 27, which lies half-way between the experimentally measured 
bounds (5 and 50). Thus, the model yields results that are consistent with the data. 

Fig. [7b is characterized by two properties: (a) There is no bistability at small ai, and (b) 
bistability is feasible at high ai, provided 6i lies lies within a finite interval. To a first degree of 
approximation, these properties reflect the following physical fact: Bistability occurs precisely when 
the induction rate is cooperative. To see this, observe that to a first degree of approximation, we 
can identify the existence of cooperative kinetics with the existence of an inflection point on the 
induction curve. Since the induction rate has an inflection point at ei = ( ai/3 — 1 1 /^i, it follows 



1. If «! is small, there is no inflection point. Hence, the induction rate is not cooperative, i.e., 
its slope decreases monotonically with ei, and bistability is impossible. 

2. If «! is large, the induction rate contains an inflection point, but it is cooperative, in effect, 
only within a finite interval of ^i. At large ^i, the inflection point is so close to that the 
induction rate is not cooperative, except at vanishingly small enzyme levels. At small 6i, the 
inflection point is so large compared to 1 that the induction rate is effectively linear on the 
interval, < ei < 1, containing the steady state. Thus, bistability is feasible only in a finite 
interval of 5i . 

In other words, the cooperativity of the quasisteady state induction rate depends on the intracellular 
parameter, ai, and the state of the environment (^i is proportional to cxi). If ai is small, bistability is 
impossible because there is no extracellular TMG level that can make the induction rate cooperative. 
If ai is large, there is a finite range of extracellular TMG levels at which the induction rate is 
cooperative, and bistability is feasible. 

3.1.2 DNA looping 

We assumed above that the repression was entirely due to repressor-operator binding. In reality, 
the repression is dominated by DNA looping. Since the induction kinetics are qualitatively different 
in the presence of DNA looping, it is relevant to ask if the foregoing conclusions are dramatically 
altered when ai > 0. We show below that the bifurcation diagram is essentially unchanged if 
< ai < 16, but it is qualitatively different at larger values of di. 

In the presence of DNA looping, the steady states are given by the equation 



For each fixed di > 0, this equation defines the steady state surface in the ai^iei-space. The folds 
on the steady state surface also satisfy the equation 



that 



d7 = 



l + ai/(l + ^iei)^ + di/(l + 5iei) 



1 



4 - ei = 0. 



2ai6i/ (1 + ^lei)^ + 4diJi/ (l + ^lei)^ 

r n 2 

1 + ai/ (l + ^lei)^ + di/ (l + ^lei)"^ 
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It is shown in Appendix [B] that the folds have the parametric representation 

1 ai(3xi-l)-(l+Xi)' 



ei(xi) 



1 



2x1 di/(i + xir-r 

Xl _ o 2"l/(l + Xl) 



4 



5i(xi) = ^H = 2xt 



ei(xi) 1-Xi 
Furthermore, bistabihty is feasible at every ai > 0, but there are three types of bifurcation diagrams: 

1. If the repression due to DNA looping is small (ai < (5/3)^ « 12.8), the bifurcation diagram 
is similar to the one obtained in the absence of DNA looping: Bistability is feasible only if 
«! is sufficiently large; moreover, even if this condition is satisfied, bistability occurs within a 
finite interval of 5i (Fig. [T3k). 

2. At intermediate levels of repression ((5/3)^ < di < 16), the bifurcation diagram contains two 
distinct bistable regions (Fig.fTSb). However, one of these regions is so small that it is unlikely 
to be observed in practice. 

3. If the repression due to DNA looping is large (di > 16), the bifurcation diagram is qualitatively 
different: Bistability is feasible at every ai > (Figs. [8h,b). This reflects the fact that if the 
repression due to DNA looping is sufficiently large, the induction rate is cooperative even if 
there is no repression due to repressor-operator binding. 

Thus, for all practical purposes, there are only two types of types of bifurcation diagrams. 

The simulations are consistent with the data shown in Fig. [2^. In wild-type cells, bistability 
occurs over the range 10 < (5i < 130 (Fig. Eb), which is in reasonable agreement with the 10-fold 
range observed in the experiments (3-30 in Fig. [2^). 

At parameter values corresponding to the cells transfected with the lac operator (ai 0, 
di 30), there is bistability, but the range of extracellular concentrations supporting bistability 
is so small that the system is practically monostable (Fig. [8li). However, the model cannot be 
compared to the data shown in Fig. [2b, since it does not account for the diffusive flux of TMG, 
which is significant at extracellular TMG concentrations exceeding 50 /iM. 

We can estimate Si^m by comparing the model prediction with the experimental data. Figs. El; 
and [2^ s how that at the upper fold point, 5i 140, and si^o ~ 30 /uM, respectively. Since Ks^i = 
680 JOzbudak et al.l . l2004l ) 



- - 30 
130 « (5i = (^1 m ^ m ~ 3100. 

'680 + 30 

This estimate is consistent with the experimental data. Indeed, the model implies at saturat- 
ing concentrations of T MG, the intracellular concentration of TMG is K~^Si^m- Since ~ 7 
HM (lOehler et al.l . l2006l . Fig. 4B), the intracellular TMG concentration at saturating condi tions is 



20mM, which is in reasonable agreement with the experimentally measured value of 15 mM (IKepe 



19601 . Fig. 5). 
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Fig. 8: Upper panel: Bifurcation diagrams for growth on TMG of (a) wild-type cells (ai = 1250), 
and (b) low-repression cells (di = 30). Lower panel: Variation of the steady state enzyme 
activity (black curve) and fluorescence intensity (green curve) with the extracellular TMG 
concentration (5i) for (c) wild-type cells (di = 1250, ai = 20), and (d) low-repression cells 
(di = 30, ai = 0). 



3.2 Growth on lactose 
3.2.1 No DNA looping 

In the presence of lactose, the steady states satisfy the equation 

^(^i) ^ 1^ /n^, ,2 - = 0, (35) 
1 + ai/ [l + diei) 

which defines the steady state surface in the ai(5iei-space. 

It is clear that no matter what the parameter values, there is at least one steady state. Indeed, 
([35]) implies that the net rate of enzyme synthesis, ^(ei), is positive if ei = 0, and negative if ei is 
sufficiently large. It follows that there is at least one 1/(1 -|- ai) < ei < 1 at which g{ei) is zero. 

It turns out, however, there is exactly one steady state because multiple steady states are 
impossible. To see this, observe that multiple steady states are feasible only if the steady state 
surface contains singular points, i.e., there are points ai, 5i, ei > satisfying ((35]) and the necessary 
condition 

9eM) = - ^^-2ei=0. (36) 



l + ai/{l + 5ieif 



3 Results and Discussion 



18 




^1 ^1 

Fig. 9: Intuitive explanation for the absence of bistability during growth on lactose of cells lacking 
DNA looping. The enzyme synthesis and dilution rates are represented by full and dashed 
lines, respectively. Stable and unstable steady states are represented by full and open circles, 
respectively, (a) During growth on TMG, bistability is feasible because the dilution rate is 
proportional to the enzyme level, (b) During growth on lactose, bistability is infeasible 
because the stabilizing effect of dilution (ef ) is so strong that an unstable steady state (and 
hence, bistability) is infeasible. 



But there are no such points because g^^ < at every point on the steady state surface. Indeed, 
([35]) implies that every point on the steady state surface satisfies the relations 



l + ai(l + 5iei) 



Substituting these relations in ([361) yields 



1 



1. 



ffei(ei) = 2ei 



1 



6iei 



1 



< 0. 



(37) 



It follows that there are no singular points on the steady state manifold, and hence, no multiple 
steady states. 

We pause for a moment to give the physical meaning of the foregoing mathematical argument. 
To this end, observe that bistability can occur only if the system permits the existence of an unstable 
steady state, i.e., a steady state at which the slope of the induction rate exceeds the slope of the 
dilution rate (Fig. [3i)Il Such a steady state is feasible when the cells are grown in the presence 
of TMG because the dilution rate increases linearly with the enzyme level. It is infeasible during 
growth on lactose because the dilution rate, ef, increases so rapidly with the enzyme level that at 
every conceivable steady state, the slope of the dilution rate exceeds the slope of the induction rate, 
i.e., < at every possible steady state (Fig. [9b)- In other words, the stabilizing effect of dilution 
is so strong that an unstable stable steady state, and hence, bistability, is impossible. 

The model assumes that the yield is constant. It is well known, however, that the yield is 
vanishingly small un der starvation conditi ons, and increases progressively as the nutritional status 



of the cells improves JTempest et a,l.l . ll967l V Thus, it seems plausible to assume that at low enzyme 



^ The stability of a steady state is completely determined by the relative slopes of the induction and dilution rates 
at the steady state. If the slope of the induction rate at a steady state exceeds the corresponding slope of the dilution 
rate (see open circle in Fig.[9ti), the slightest increase (resp., decrease) in the enzyme level increases (resp., decreases) 
the induction rate more than the dilution rate, so that the enzyme level grows (resp., declines) even further. 
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Fig. 10: Growth of wild-type cells on lactose: (a) Bifurcation diagram (di = 1250). The dashed line 
shows the path along which 6i changes when the extracellular lactose concentration is varied, 
(b) The variation of the steady state enzyme activity (black curve) and fluorescence intensity 
(green curve) with 5i, a surrogate for the extracellular lactose concentration (di = 1250, 
ai = 20). 



levels, the lactose transport rate and yield are relatively small; as the enzyme level increases, so 
do the lactose transport rate and the yield. If this is true, it is relevant to ask if multistability 
is impossible even if the yield is an increasing function of the enzyme level, i.e., the yield has the 
form Yi(;/)(ei), where (/'(ei) < 1 is some increasing function of ei. The foregoing physical explanation 
suggests that multistability is, a fortiori, impossible because the dilution rate now increases with ei 
at a rate even faster than . Analysis confirms this intuitive argument — the value of g^^ at any 
point on the steady state surface is 

-0,,(ei)e? <0. 

Thus, multistability is impossible even if the yield is an increasing function of the enzyme level. 

The above argument shows that the existence of bistability is not determined solely by the 
intensity of the destabilizing positive feedback generated by induction. It also depends on the 
strength the stabilizing effect of dilution. If this stabilizing effect is sufficiently large, it can neutralize 
the destabilizing effect of positive feedback. 

3.2.2 DNA looping 

We have shown above that in the absence of DNA looping, there is no bistability during growth 
on lactose because the dilution rate, which is proportional to ef , neutralizes the destabilizing effect 
of positive feedback. However, in the presence of DNA looping, the destabilizing effect of positive 
feedback is much stronger because the induction rate contains terms proportional to ef. Under these 
conditions, it seems plausible to expect that bistability is feasible at sufficiently large di. We show 
below that is indeed the case, but the bistable region is extremely small. 



20(ei)ei 



{l-<^(6l)6?}- 



5iei 
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In the presence of DNA looping, the bifurcation points satisfy the equations 

1 2 _ 

^ l + ai/{l + 5ieif + ai/{l + 5ieit " " °' 

9eAei) = -2 - 2ei = 0, 

1 + Qi/ (1 + Sieif + di/ (1 + (5iei)^ 

For each fixed di > 0, these two equations define the locus of the bifurcation points in ai^iei-space. 
It is shown in Appendix O that this curve has the parametric representation 



1 



^Xi{ai/(l + Xi)'-l} 



(1+Xi 



i2 



. . X Xl 3/2 Oil 

Furthermore, bistability is feasible only if di > 5^/2^ « 195. 

In wild- type cells, for instance, bistability occurs for all 7 < 5i < 9 (Fig. [TOh). Now, it is 
conceivable that bistability is not observed during growth on lactose because the strength of the 
positive feedback, is so small (<7) that 5i = remains below the bistable region at all 

concentrations of extracellular lactose. It is clear, however, that even if 6i^rn is large (>9), the 
bistable region is so narrow that it is unlikely to be observed in practice (Fig. [TOb). We show below 
that the width of the bistable region increases dramatically in the presence of glucose. 

In wild-type cells, bistability is feasible during growth on lactose, but the bifurcation diagram is 
qualitatively different from that obtained during growth on TMG. Indeed, during growth on TMG, 
bistability is feasible for all ai > (Fig.Eb). In sharp contrast, bistability is feasible during growth 
on lactose only if ai is sufficiently small (Fig. [TOb). The qualitatively different bifurcation diagram 
for lactose can be explained as follows. During growth on lactose, bistability can occur only if DNA 
looping has a strong effect on the induction rate. But DNA looping dominates the induction rate 
precisely if 

di ai / di 

»7: ^^xi«a/ 1- 



(l + Xi)' {l + Xif V ai 

It follows that regardless of the value of di, the range of inducer concentrations at which DNA 
looping is dominant vanishes at a sufficiently high value of ai. Consequently, the induction rate, 
and hence, the dynamics, become similar to those observed in the absence of DNA looping. 

3.3 Growth on lactose + glucose 

In the presence of glucose and lactose, the steady states satisfy the equations 

gi{ei,e2) = — i^-/n^A .4 " (^i + "^2) ei = 0, (38) 

I + ai/ [1 + Oiei) +ai/{l + diei) 

o 

52(ei,e2) = -— 777— T — r - (ei + ae2)e2 = 0, (39) 

1 + "2/ (1 + 0262) 
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Fig. 11: Classification of the dynamics during growth on lactose + glucose. Bistability is feasible 
if and only if ai and ai lie above the full curve. If ai and ai lie between the full and 
dashed curves (resp., aboe the dashed curve), the bifurcation diagram has the form shown 
in Fig. m (resp., Fig. [I2ji). 



Our goal is to address the following question: Given a particular cell type, what are concentrations 
of lactose and glucose at which bistability is feasible? This question is difficult to address because 
there are three parameters that depend on the substrate concentrations {6i oc a oc \/o^Jo\). 
Since the combination of parameters 



(5 = 



a6i 



-1 

x,2 



xA 



is independent of the substrate concentrations, it is convenient to replace a with f362 /6i. In terms of 
the model, the question of interest then becomes: Given any fixed ai,a2, Ci2, P, what are the values 
of Si and S2 at which bistability is feasible? Unlike 01,02,02, the value of /? cannot be determined 
from the experimental literature, since K~2 is a phenomenological parameter. In the simulations, we 
assume that /3 = 1. However, we show below that this is not a particularly restrictive assumption, 
since the qualitative behavior of the bifurcation diagram is completely determined by ai and ai. 
Bistability is feasible only if there are steady states satisfying ([38]) - ([39] l and the condition 



det 



dgi dgi 

dei de2 

§92 dg2 

dei dt2 



It is shown in Appendix [D] that these three equations can be rewritten as 



h{xi,X2) = Xi V 1 + "1/ (1 + Xi) + ai/ (1 + Xi) 



S2{Xl,X2) = X2\/l + a2/(l + X2)l 



P = —Hxi,X2), 

X2 



1 + ^(Xl,X2) 
^(Xl,X2) 



1 + ^(Xl,X2) 



(40) 

(41) 
(42) 
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(a) (b) 

Fig. 12: Steady states of wild-type cells during growth on glucose and lactose [ai = 20, ai = 1250, 
a2 = 4, /3 = 1). (a) Bifurcation diagram. Bistability occurs if 6i, 62 lie in the region between 
the curves, (b) Variation of the steady state activities of Ei (black curve), E2 (red curve) 
and G (green curve) with the extracellular lactose concentration, 61, at the fixed glucose 
concentration, S2 = 10. 



where 

l/p(Xi) - 1 

1 - q{x2) 

^ ai/(l + Xi)' + 2ai/(l + Xi)^ Xi ^ 
l + ai/(l + Xi)' + ai/(l + Xi)'l + Xi ' 

Q2/(l + X2) X2 _ ^ 

l + a2/(l + X2) 1 + X2 

Furthermore, bistability is feasible if and only if ai and ai lie above the full curve shown in Fig. [TTl 
It is precisely for such values of ai and ai that eq. (l42l) has positive solutions, which determine a 
curve in the xiX2-space. The variation of 5i(xi)X2) and (^2(xi)X2) along this curve circumscribes 
the bistability region on the (^i52-plane (for the cell type defined by the fixed values of ai,ai,a2, (3). 
If ai and ai lie between the full and dashed curves of Fig. [TTl the bifurcation curve has the form 
shown in Fig. [14] — it does not intersect the 5i-axis. If ai and di lie above the dashed curve of 
Fig- [m the bifurcation curves intersect the (^i-axis (Fig. [T2k). Thus, the form of the bifurcation 
curve is completely determined by ai and ai. 

Fig. [TT] implies that in cells lacking DNA looping, bistability is feasible if and only if ai exceeds 
the threshold value of 27, which is identical to the threshold in such cells when they grow in the 
presence of TMG (Fig. [7|. The existence of this identity is not coincidental. As we show below, the 
dynamics of growth on lactose + glucose are, in some sense, identical to the dynamics of growth in 
the presence of TMG. 

Fig. [T2k shows the bifurcation diagram for wild-type cells growing on lactose + glucose. Evi- 
dently, bistability is feasible at every 62 > 0. This mirrors the fact that during growth of wild-type 
cells on pure lactose {62 = 0), bistability is feasible, although the range of lactose concentrations 
supporting bistability is very small. At relatively modest values of the extracellular glucose concen- 
tration, 62, the bistable region increases dramatically. This is because at sufficiently large concen- 
trations of glucose, the dilution rate due to growth on glucose becomes significant. Importantly, this 
dilution rate is formally similar to the dilution rate during growth on TMG, insofar as it increases 



^(Xi,X2) = 

p{xi) = 
q{x2) = 
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linearly with ei. Thus, the dynamics are also analogous to those observed during growth on TMG: 
Bistability is feasible for a range of lactose concentrations {6i). 



Fig.[T2b is qualitatively similar to the bifurcation diagram obtained by Santillan et al (ISantillan et al 



20071 . Fig. 2). However, in their model, the width of the bistability region increases in the presence 



of glucose due to regulatory effects, such as cAMP activation and inducer exclusion, exerted on the 
induction rate of the lac operon. In our model, the same phenomenon occurs because the stabilizing 
effect of dilution is significantly smaller during growth on lactose + glucose (as opposed to growth 
on pure lactose). 

It remains to explain the discrepancy between the data obtained Ozbudak et al and Loomis & 
Magasanik. The simplest explanation is that in the strain used by Ozbudak et al, the parameter, 
^i,m, which measures the strength of the positive feedback, is so small (<7) that 5i = (^i,m.ci does 
not enter the bistability region at any extracellular lactose concentration (Fig. [T2b). On the other 
hand, the values of 5i^m, ^2,m for the strain used by Loomis & Magasanik are so large that 81,82 lie in 
the bistable region at the extracellular lactose and glucose concentrations used in their experiments. 

It should be noted, however, that the discrepancy could also reflect technical differences in the 
experiments. Fig [T2b shows the steady state levels of the lactose enzymes (black curve), glucose 
enzymes (red curve), and GFP (green curve). In the neighborhood of the upper bifurcation point 
{81 « 35), the lactose enzyme level of the induced cells is ~100 times the lactose enzyme levels of 
the non-induced cells. However, the GFP level of the induced cells is no more than ~ 10-20 fold 
higher than the GFP level of the non-induced cells. Now, Fig. [2^ shows that in the neighborhood of 
the upper bifurcation point, the steady state fluorescence of the non-induced cells is scattered over a 
~50-fold range (see the fluorescence distribution in lower panel of Fig. [2^ over the range 20-30 /xM 
). Assuming that the distribution has a similar variance during growth on lactose, it is conceivable 
that even if bistability exists, the fluorescence distributions of the non-induced and induced cells 
overlap, and appear to be unimodal. 

In Fig [T2b . the precise values of the lac enzyme and GFP levels in induced and non-induced 
cells depend on the particular choice of the parameter values. However, as shown above, the ratio of 
GFP levels in induced and non-induced cells is always smaller than the corresponding ratio for the 
lac enzymes. Experimental artefacts indicating the absence of bistability are therefore more likely 
if GFP levels, rather than enzyme activities, are measured. 



4 Conclusions 



The experimental data shows that the existence of bistability in the lac operon depends on the 
composition of nutrient medium. It occurs during growth on TMG/succinate, but not on lactose. 
There are conflicting reports of its existence in media containing lactose and glucose. 

The occurrence or absence of bistability reflects the net result of the destabilizing effect due 
to the positive feedback generated by induction and the stabilizing effect of dilution. In previous 
models, the experimental data has been rationalized entirely in terms of changes in the strength of 
positive feedback as a function of the medium composition. We have shown above that: 

1. The stabilizing effect of dilution also changes dramatically with the composition of the medium. 
In the presence of non-galactosidic carbon sources, such as succinate or glucose, the dilution 
rate of the lac enzymes contains a term that is proportional to the activity of these enzymes. 
During growth on pure lactose, the dilution rate of lac enzymes is proportional to the square 
of their activity. 

2. These variations in the functional form of the dilution rate have a profound effect on the dy- 
namics. During growth on lactose, the stabilizing effect of dilution is so strong that bistability 
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is virtually impossible even if induction is subject to the strong positive feedback generated by 
DNA looping. During growth on TMG/succinate or lactose + glucose, bistability is feasible 
because the stabilizing effect of dilution decreases sharply. Thus, bistabililty is much more 
likely in the presence of non-galactosidic carbon sources, namely, succinate and glucose. 

3. The conflicting results on bistability during growth on lactose + glucose can be explained in 
terms of the relative magnitudes of the destabilizing and stabilizing effects of positive feedback 
and dilution, respectively. However, the criterion used by Ozbudak et al to infer monostability, 
namely, the absence of a bimodal distribution in non-induced cells, is prone to error. Since 
the repression of the reporter lac operon is much lower than the repression of the native lac 
operon, the ratio of GFP levels in induced and non-induced cells is significantly lower than 
the corresponding ratio of the enzyme activities. The fluorescence distribution can therefore 
appear to be unimodal even if the enzyme levels are bistable. 

Taken together, these results show that while the intensity of the positive feedback undoubtedly 
influences the dynamics of the lac operon, the dilution rate also has profound effects. These effects 
can be discerned only if it is recognized that the specific growth rate is not necessarily a fixed 
parameter — it depends on the physiological state of the cells. 
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A Derivation of eqs. (fT8])-([2T]) 

Since xi + X2 + ei + g + 62 + C = 1, addition of eqs. (fT3]) - (fT7|l yields 
where 



Idc 
c dt ' 



E 

1=1 



Ks.i ~\~ Si 



K,^^^ - (1 - Yi) /C+ 



Xi 



is the specific growth rate. This becomes evident if we rewrite the above equation in the more 
familiar form 

dc 

di = 

Alternatively, one can see that as expected, rg is the net rate of uptake of the two substrates (uptake 
minus loss by excretion and respiration). 

If the experiments are started with extremely small inocula and terminated before the cell 
densities become significantly large, the substrate concentrations do not change significantly over 
the course of the experiment, i.e. Si{t) « Si^ = Sj(0) for all t. Since fcj'j > k~- is large, Xi rapidly 
attains quasisteady state, i.e., eq. (fT3l) becomes 



K.,:e,; 



Sjfl 

Ks.i + s. 



i,0 
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which imphes that 



Vs,i _ Sjfl 

— r ■ T^eiai, (Ji = — ■ , 

2 2 

^.^iK.i^i = '^Yi(j)iVs^ieiai, (f), 
i=i ' j=i 



(43) 
(44) 



Eq. (03]) is identical to eq. Since {l/c){dc/dt) = rg, substituting ^ in ([I4|-([T6|) yields 

m-m. 



B Bifurcation diagram for growth in the presence of TMG 

The bifurcation points satisfy the equations 

1 

ei = 0, 



fiei) 



1 + ai/ (1 + (^lei)^ + di/ (1 + 6iei)^ 
2ai6i/ (1 + Sieif + 4di^i/ (l + ^lei)^ 



1 + ai/ (1 + Sieif + di/ (1 + dieiY 



1 = 0, 



which imply that 



+ 



ai ^ Qi 



+ 2- 



1 


-1, 






1 + 



1 



5ieiJ ei 



1 



Since xi = ^i^i, we have two equations that are linear with respect to ai, di, and 1/ei. Solving 
for ai and 1/ei in terms of di and xi yields 



ai(xi) 
ei(xi) 
^i(xi) 



1 



1 1-Xi 



Mxi) ^ ai (3X1 - 1) - (1 + Xif , 



2X1 V'2(xi 
Xi 



1-Xi 

, ^2(Xl 



(l + Xl 



ei(Xi) '^'1-Xi 

These three equations provide a parametric representation of the bifurcation curve. 
In the absence of DNA looping, the parametric representation of the curve becomes 

\3 



/ \ Xl - 1 . ^ (1 +Xl) . / X 

ei(Xi) = , ai(Xi) = , ^i(XiJ 



2X1 



Xi 



1 



2xi 

XI -1' 



It follows that the bifurcation curve exists (i.e. lies in the positive octant of the ai^iei-space) for 
all xi > 1- 

To determine the existence of bistability in the presence of DNA looping, we begin by observing 
that ei(xi) > whenever 5i(xi) > 0. Hence, it suffices to confine our attention to ai(xi) and 
5i{xi)- But ai and 6i are positive on the interval < xi < 1 (resp., 1 < xi < oo) if and only if ipi 
and ip2 are positive (resp., negative). Thus, we are led to consider the signs of V'l and ^^2 on xi > 0. 

It follows from the analysis of ipi that: 
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1. If < ai < (5/3)^ then -01 < for all xi > 0. 

2. If ai > (5/3)^, then ipi has two positive roots, < < Xii and tpi > if and only if 
Xi < Xi < Xi- Furthermore, XiiXi < 1 when (5/3)^ < ai < 16, and X^^ < 1 < Xi when 
ai > 16. 

Similarly, the analysis of ■02 shows that: 

1. If < ai < 1, then V'2 < for all xi > 0. 

2. If ai > 1, then ^^2 > if and only if < xi < X* = ^h.^^ ~ ^- Furthermore, if oi > (5/3)^, 
then Xi lies between Xi and 1. 

Taken together, these results imply that there are three distinct types of the bifurcation diagrams: 

1. If < ai < (5/3)^, then ai(xi)) ^i(xi) are positive if and only if xi > 1 (Fig. [T3b). 

2. If (5/3)^ < ai < 16, then ^i(xi)) cki(Xi) are positive if and only if X;^ < Xi < Xi or x > 1, 
where < x^^ < Xi < 1 are the positive roots of ai(xi). Each of these two intervals of 
existence yields a bistability region (Fig. [T3b). 

3. If ai > 16, then 5i(xi); oi(xi) are positive if and only if x^ < X < 1 or x > Xi, where Xi and 
Xi are the positive roots of ai(xi) (Fig. [8h,b). 

The geometry of the cusps in the foregoing figures follow from the relations 

dai _ ^ ai {3x1 - 3X1 + 2) + (1 + xi)' (xi - 2) 

dXi (l+xi)'(l-Xi)' 
d5i xi dai 

dxi {l + Xifdxi' 

which imply that there is a cusp whenever dai/dxi = 0. It follows from the properties of ■0i that 
there are two cusps if and only if (5/3)^ < ai < 16. If ai > 16, one of these cusps disappears 
because its cti-coordinate is negative. 
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C Bifurcation diagram for growth on lactose 

In this case, the bifurcation points satisfy the equations 



I ,2 

l + ai/(l + <5iei)2 + ai/(l + (5iei)^ ^ 
2ai6i/ (1 + 5ieif + 4ai5i/ (1 + 6ieif 



1 + ai/ (1 + 6ieif + ai/ (1 + 6iei)^ 



2ei = 0. 



which imply that 
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(l + ,5iei)^ (l + <5ieif 
ai ai 



+ 2- 
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1 \ 1 



Since xi = ^i^i , we have two equations that are hnear with respect to ai, cii, and 1/ef. Solving 
for ai and 1/ei in terms of ai and xi yields 



ei(xi) 

aiixi) 
Siixi) 



^ Xi{«i/(l + Xi)^-l}' 



e(xi 



2 ' 



V'3(xi) = "1 ixi - 1) - (1 + xi)" 




Ql 



For each fixed ai > 0, these relations provide a parametric representation of the bifurcation curve. 

In the absence of DNA looping, the bifurcation curve does not exist because the ei- and 6i- 
coordinates of the curve are imaginary (and the ai-coordinate is negative) for all xi- 

In the presence of DNA looping, the bifurcation curve exists for all ai > 5^/2^ 195. To see 
this, observe that since 



ai(xi) > ^ ai (xi + 1) > ai (xi - 1) > (1 + Xi)' 



(45) 



ei and 5i are positive whenever ai is positive. Hence, it suffices to confine our attention to ai. Now, 
ai > if and only •i/'s > 0. One can solve the equations ■03 = dips/dxi = to conclude that ips > 
for some xi > if and only if di > 5^/2^. In this case, ip3 and ai have two roots, < x^ < Xi) 
and are positive if and only if xi lies between these roots. 

The bifurcation diagram is qualitatively similar to Fig. [TOk for all di > 5^/2^ because under 
this condition, ai and 6i are positive on {x^,Xi)j and achieve a unique maximum on (x^^^Xi) at the 
very same value of xi- To see this, observe that ai has at least one maximum in (x^^jXi)- In fact, 
it has exactly one maximum because 

doi _ ai (3 - xi) - 3 (1 + xi)^ 
dXi~ {l + Xif 
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cannot have more than one zero on this interval. On the other hand, 5i > on [x^^iXi] because 
6i > «!. Furthermore, 5i and ai attain a maximum at the very same Xi because 

d6i 1 / xi \^ 1 '^Q^i 
dxi 2V1 + Xiy Sidxi' 
(f6i 1 (Pai 

The first relation implies that d6i/dxi = precisely when dai/dxi = 0. The second relation implies 
that when d5i/dxi (and hence, dai/dxi) is zero, d'^di/dxf and d'^ai/dxf < have the same sign. 

The width of the cusp-shaped region always increases as one moves away from the cusp point 
(Fig. [TOk). To see this, observe that the above equations imply that 



dai (1 + xi)^ 

is an increasing function of xi- Hence, the slope at any point on the upper branch of the bifurcation 
curve is always greater than the slope of any point on the lower branch. 



D Bifurcation diagram for growth on lactose + glucose 

In this case, the steady states satisfy the equations 

Pi = (ei + ae2) ei, pi 



1 



P2 = (ei + ae2) e2,P2 



1 + ai/ (1 + dieiY + ai/ (1 + 6iei) 
a 



l + a2/(l + (52e2)' 



(46) 
(47) 



where a = /35i/52- If a steady state corresponds to a fold bifurcation point, the determinant of the 
Jacobian, 



J 



2ai<5i/(l+,;iei)='+4ai<5i/(l+<5iei)'^ „^ 

T o TTT — '^^1 — 0:t2 



-£2 



~aei 

a 0262 _ , _ 2fyf 

{l+a2/(l+<52e2)}^ {1+S2e2f ^ ^ 



at that steady state must also be zero. Now, it follows from ([46]) that at a steady state 

Jn = (ei + ae2 ei— — — — — - 2ei - ae2 

1 + ai/ (1 + <5iei)' + ai/ (1 + <5iei)^ 

= (ei + 062) p- ei 



where 



P 



6i€i ai/il + 6ieiy + 2ai/{l + 6iei) 



'1 + Siei 1 + ai/ (1 + 6ieif + di/ (1 + 6iei)^ 
Similarly, ([47|l implies that at a steady state 

1 2 



1. 



J22 = (ei + 0^2) £2— , .1 , . :r- 72 

1 + 02/ (1 + S2€2) (1 + 6262) 

= (ei + 062) q - ae2, 



ei - 2a€2 
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where 



a2 



h^2 



"2 + (1 + (5262) 1 + ^2^2 

It follows that det J is zero at a steady state if and only if 



1. 



I.e., 



(ei + 062) [(ei + 062) P9 - a-pe2 - ?ei] = 0, 

a = —h, h = — . 

£2 1 - 1/q 



(48) 

The bifurcation points satisfy eqs. (pSil - ipHI) . 

To determine the parametric representation of the bifurcation points, observe that (l46]l - (l47|l 
yield ei/e2 = P1/P2, which can be substituted in ([48]l to obtain a = {pi/ P2) h, i.e., 



a 



1 + Q2/(1+X2) 



where h is now a function of xi and X2- Eqs. (l46]) - (l47ll also imply that 



Vh, 
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Finally, 
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It follows from ^M that 



£2 



Vl + Q2/(1 + X2) V l + h' 



and 



'52(Xl,X2) = — = X2\/l + 02/ (1 + X2) 
£2 
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Substituing (I50l)-(l5l|) in ^ yields 



X2 a2 + (l + X2) P(xi) 



(49) 



(50) 
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(52) 
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Fig. 14: The bifurcation diagram for growth on glucose and lactose at ai = 170, ai = 0, 02 = 4, 
(3 = 1). The (^2-coordinate of the bifurcation curve is always positive, i.e., bistability is 
infeasible at sufficiently small 82. This reflects the fact that if the repression due to DNA 
looping is small, bistability is impossible during growth on pure lactose {82 = 0). 



Eqs. ([50]l - ([52]) provide a parametric representation of the bifurcation points. 

Bistability is feasible if and only if and only if < p < 1 for some xi > 0- To see this, observe 
that the LHS of ([52]) is a monotonically increasing function of X2 for all X2 > 0, and the RHS of 
([52]) is positive if and only if < p < 1. It follows that for every xi > such that < p{xi) < 1, 
there is a unique X2 > satisfying ([52]) . These pairs, Xi)X2 > 0, define a curve on the xiX2-plane, 
and the variation of 6i{xi,X2) and 52(xi)X2) along this curve delineates the bistable region on the 
5i52-plane. In contrast, for every xi > such that p{xi) < or p{xi) > 1, the RHS of ([52]) is 
negative. Hence, there is no X2 > satisfying ([52]) . and no bistability. 

The condition, p > for some xi > 0, is satisfied if and only if ai and di lie above the full 
curve shown in Fig. [TTl To see this, observe that p has a maximum on xi > since p{0) = — 1, 
p{oo) = —1, and p{xi) > —1 for < xi < 00. The value of ai and di at which a maximum of p 
touches the xi-axis satisfies the equations 

01 di 
Pixi) = ^ (xi - 1) ,2 + (3X1 - 1) ,2 = 1 + XI, 

Qi di 
PxAxi) = 0^{3-xi) 7-^-2 + (7 - 9xi) ,^ , ,2 =1 + ^1 

which can be solved to obtain 

Jl + Xi)'(3xi-2) . . , (l + xi)'(2-xi) 



3x1-3x1 + 2 3;^^_3^, + 2 

The above relations define the full curve shown in Fig. [TTl Evidently, the di- and ai-intercepts of 
the curve are 27 and (5/3)^, respectively. Moreover, since dp/dai > for all xi > 0, the condition, 
< p for some xi > 0, is satisfied precisely when ai and di lie above this curve. 

If ai and di lie above the full curve in Fig. [TTJ there is an interval, say, (x^^jXi)) on which 
p > 0. But there are two possibilities. Either p never exceeds 1 on this interval, or it exceeds 1 
on some subinterval of (Xj^iXi); say, (xi,/)Xi,m)- In the first case, the X2-coordinate of the curve 
defined by ([52]) . and hence, the 52-coordinate of the bifurcation curve is always positive, i.e., the 
bifurcation diagram has the form shown in Fig. [121 In the second case, the X2-coordinate of the 
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curve defined by (|52p is zero at xi — Xi,h Xi,u^ 

and the 52-coordinate of the bifurcation curve is zero 
at the corresponding point. The bifurcation diagram therefore has the form shown in Fig. I12h. 

One can check (by the method similar to the one shown above) that p exceeds 1 for some xi > 
if and only ai and Qi lie above the curve defined by the relations 

^ (l + Xi)'(2xi-3) . Jl + Xif 

aiixi) = 2 , ai(xi) = 3^ , 

3 - Xl 3 - Xi 

which is shown as the dashed curve in Fig. [TTl Evidently, the di- intercept of the curve is 5^/2^. 
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